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SUMMARY 


A fluid mechanical model of the acoustic behavior of small 
orifices is presented which predicts orifice resistance and 
reaotanoe as a function of incident sound pressure level, frequency, 
and orifice geometry. Agreement between predicted and measured values 
(in both water and air) of orifice impedance is excellent. The 
model shows that 

(1) The acoustic flow in the immediate neighborhood of the orifice 
(i.e., the near field) can be modeled as a locally spherical flow. 
Within this near field, the flow is, to a first approximation, 
unsteady and incompressible. 

(2) At very low sound pressure levels, the orifice viscous resistance 
is directly related to the effect of boundary -layer displacement 
along the walls containing the orifice and the orifice reactance is 
directly related to the inertia of the oscillating flow in the ori- 
fice neighborhood. Previously, orifice resistance and reactance 
were modeled by empirical end correction expressions. The model 
also shows that at low to moderate sound pressure levels, the resis- 
tance can be dominated by weak nonlinear jet-like losses but that 
the overall impedance can still be constant (i.e., independent of 
incident sound pressure level) providing the orifice resistance is 
very small relative to the reactance. This is shown to occur when 
the amplitude of the incident acoustic pressure P is less than 

p [to (D+L) ] ^ where to is the sound radian frequency, (D+L) is the ori- 
fice diameter and thickness respectively and p is the fluid mean 
density. 

(3) When P/p [to (D+L) ] ^>>1, the orifice impedance is dominated by 
nonlinear jet-like effects. This corresponds to very high sound 
pressure levels at which the orifice behaves in a predominately 
quasi- steady manner. Thus the model establishes explicitly the 
quasi-steady nature of the flow in orifices exposed to intense sound. 

(4) When P/p [to(D+L) ] ^=0 (1) , orifice resistance and reactance are 
roughly equal. 
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DEFINITION OF SYMBOLS 
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Definition 

constant defined by Eqn. (36) 
airy function 
speed of sound 

time -averaged discharge coefficient 
orifice diameter 

special functions defined in text 
special function defined in text 
orifice thickness 

radial distance of truncated hemispherical 
surface corresponding to orifice outlet 

acoustic pressure 

amplitude of incident acoustic pressure 
radial distance 
acoustic resistance 

orifice Reynolds number (o)* (D* + L’^) v*) 

time 

radial velocity 

maximum velocity at orifice vena contracta 
orifice inlet velocity 
orifice reactance 
orifice impedance 

parameter defined as (D*+L*) /V* 

parameter defined as V*/o3* (D*+L*) 

parameter defined as /2~3 

boundary- layer displacement thickness 

ratio of specific heats 

fluid kinematic viscosity 

fluid coefficient of viscosity 

fluid density 

sound wavelength 

sound radian frequency 

spherical polar angle (See Fig. 4a) 

angle defined as tt/2-0 

transformed time coordinated (Eqn. 99) 


V 



Symbol 


Definition 



Superscripts 

C )* 

n* 


transformed boundary layer coordinate (Eqn. 54) 
nonlinear orifice end correction (Eqn. 2) 


denotes dimensional quantities 
denotes mean quantities 


Subscripts 

0 denotes lowest-order term 

1 denotes first-order term; also harmonic term 

where obvious 
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INTRODUCTION 


Cavity backed orifices are extensively used in the aircraft 
industry as acoustic devices to reduce or absorb internally generated 
jet engine machinery noise. The efficient application of these de- 
vices depends intimately upon the selection of the ’’optimum*' impe- 
dance to maximize the sound absorption. The sound absorption theories 
of Morse ^ and Cremer^ for rectangular ducts without flow show that 
the sound absorption decreases rapidly from its maximum value for 
off -optimum wall impedance. This sensitivity has also been shown by 
Rice^ to exist for ducts containing flow. These studies demonstrate 
clearly the importance of accurately specifying the wall impedance 
in acoustically treated ducts. 

Despite the extensive use of cavity-backed orifices in industry 
as devices to absorb undesired sound, their detailed acoustic behavior 
is not well understood. It has been shown by Ingard^ and others that 
the absorption characteristics of these devices are directly related 
to their impedance. Thus, most acoustic studies of the behavior of 
cavity-backed orifices consist of the measurement and prediction of 
their impedance. The purpose of this report is to present a fluid 
mechanical model of the behavior of isolated small orifices as a 
function of incident sound pressure level, frequency, and orifice 
geometry. It is believed that this model will provide the necessary 
first step in understanding the behavior of cavity-backed orifices. 

Rayleigh® was the first to predict the impedance of orifices 
by using the concept of lumped elements in a simple mechanical os- 
cillator analogy (i.e., the slug-mass model). His model is essentially 
non- fluid mechanical but gives the actual acoustical impedance char- 
acteristics for low sound pressure levels when an empirical end cor- 
rection is added to the slug mass. Rayleigh’s model was modified 
first by Sirignano® and later by Zinn^ by introducing fluid mechanical 
concepts. To simplify their models, they assumed that the character- 
istic dimensionsof both the orifice or cavity are very much smaller 
than the incident acoustic wavelength and, further, that the acoustic 
flow through the orifice is one-dimensional, incompressible , quasi- 
steady, and calorically perfect. 

Both authors base their models on an integral formulation of 
the conservation of mass and momentum applied to two control volumes, 
one being the volume bounded by the orifice inlet and outlet surfaces 
and the other the cavity. To solve these integrals, they used the 
method of successive approximations with the first order solution 
corresponding to the linear case of very small sound pressures inci- 
dent to an orifice. The orifice nonlinear behavior is introduced 
through the higher order terms and represent only a second order 
approximation to the (linear) first order solution. Thus their 
conclusions apply only to weakly nonlinear acoustic pressures and not 
to the intense sound pressures existing within rocket chambers or 
jet engines, the intended application of their models. 

There is a serious deficiency common to both of their models. 
Sirignano assumes the loss in acoustic energy at the orifice outlet 
is equal to the jet outlet kinetic energy. Zinn assumes that, at 



the orifice inlet, the axial inlet flow is zero but allows a radial 
inflow to preserve continuity. Both of these assumptions are 
difficult to understand because they violate their original assump- 
tions. For example, Zinn's assumption that the flow in the orifice 
is one-dimensional (i.e., 9u/9x = 0) clearly contradicts his zero 
inflow and jet-like outflow assumption. Sirignano violates the 
conservation of momentum by arbitrarily including a momentum term 
equal to ypu^/2 (see the third term on the RHS of Eqns. 10(a) and 
10(b) of his paper). It is interesting to note that these two assump 
tions lead directly to a one-half difference in their estimate of the 
orifice nonlinear resistance. Another major deficiency of their 
models is that to first order (i.e., the so-called linear orifice 
impedance regime), both models predict the cavity resistance but not 
the orifice reactance. 

Despite these criticisms, Sirignano and Zinn were .the first to 
assume that the behavior of the nonlinear acoustic flow in the 
neighborhood of the orifice is quasi-steady and that the concept 
of a discharge coefficient properly connects the orifice inflow to 
the outflow. 

Measurements of the behavior of small isolated orifices by 
Ingard and Ising® and by Thurston® et . al have provided valuable data 
and much needed physical insight. These studies are reviewed below 
because of their importance in the development of the fluid mechani- 
cal model described in Section 2. 

Ingard and Ising used the arrangement shown in Figure 1 to 
study experimentally the acoustic nonlinearity of an isolated orifice 
Figure 1 shows an orifice plate mounted at one end of a circular 
cylinder. The experimental program consisted of taking simultaneous 
measurements of the acoustic pressure within the cavity and the 
acoustic velocity in the orifice. The sound pressure level within 
the cavity was measured with a small condenser microphone. The 
acoustic velocity in the orifice was measured with a hot-wire probe 
placed at the center of the orifice. The air within the cavity was 
excited at a frequency of 150 Hz by means of an electromagnetically 
driven piston located at the bottom of the cavity. 

The acoustic nonlinearity is described in terms of the behavior 
of the orifice impedance. Ingard and Ising defined the orifice 
impedance as the ratio of pressure within the cavity to the funda- 
mental harmonic component of the orifice inlet velocity. The funda- 
mental harmonic component of the orifice velocity was calculated by 
performing a Fourier decomposition of the measured orifice velocity 
time-history. The magnitude of the impedance is given by 

( 1 ) 

where p? and ut represent the amplitudes of the harmonic cavity 
pressure and orifice velocity respectively. The phase angle c|)i 
between the acoustic pressure p\ and velocity Ui was determined 
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graphically by comparison o£ the pressure and velocity traces dis- 
played simultaneously on an oscilloscope. 

The results o£ their study are summarized in Figure 2 in terms 
o£ R; the ori£ice resistance, and X*' the ori£ice reactance. The data 
plotted against ori£^ice inlet velocity, m^ be divided into two 
regions, one where and the other f5>>X*. In the region where 

the ori£ice impedance ^=*)C^is, to £irst order, constant indepen- 
dent o£ incident sound pressure level (and hence independent o£ 
ori£ice inlet velocity) . In the region where and the data 

shows that, to £irst order, 

Ingard and Ising o££er the £ollowing interpretation o£ their 
data. At low sound pressures (corresponding to R^<x) , the ori£ice 
resistance and reactance are given by the £ollowing empirical expres- 
sions 


where L* is the orifice plate thickness, if the diameter, the 

nonlinear resistive end correction, and 0.85D*is the reacifive mass 
end correction. These equations predict quite well the orifice 
impedance at these low sound pressure levels. The orifice inflow 
and outflow is essentially irrotational ; the acoustic driving pres- 
sure is balanced primarily by the instantaneous (i.e,, local) accel- 
eration of the acoustic velocity in the vicinity of the orifice. 

At high sound pressures (where , the measurements show 

(see Fig, 2) that the orifice resistance is proportional to the ori- 
fice velocity. The measurements also showed that the orifice reac- 
tance is very insensitive to orifice velocity, decreasing at the 
very highest sound pressure levels measured, to a value roughly one- 
half the linear value. Ingard and Ising interpreted the orifice 
resistance data in terms of Bernouilli’s Law suggesting that the flow 
behavior through the orifice is quasi- steady . The hot-wire measure- 
ments indicated that at these high sound pressure levels, the flow 
separates at the orifice forming a high velocity jet. Thus during 
one-half cycle, the flow incident to the jet is irrotational; it is 
highly rotational (in the form of jetting) after exiting from the 
orifice. During the other half of the cycle, the flow pattern is 
reversed. The loss of one-half of the reactance at these high pres- 
sure levels was accounted for by assuming that one-half of the end 
correction is ’’blown" away by the exiting jet (in their experiments 
LfD?<l hence from Eqn. 2 most of the reactance is due to the end 
correction). Ingard and Ising also measured the orifice velocity 
as a function of axial distance from the orifice (and away from the 
cavity - see Fig. 1). They found that the inflow velocity rapidly 
decayed to very small values at distances of about two to three dia- 
meters indicating that an acoustic near field existed. 
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In an earlier study, Thurston et. al measured the impedance 
of an isolated square-edged orifice of diameter 0.305 cm and, 
thickness 0.0126 cm immersed in water at temperature 26®C. The 
experimental set-up is quite similar to that used by Ingard and 
Ising consisting of a slender cylinder with a piston driver located 
at one end and the orifice at the opposite end. The acoustic 
pressure in the cavity was monitored with a capacitive type pressure 
detector and the acoustic velocity, generated by a piston-driver, 
was measured by means of a calibrated velocity pick-up attached to 
the drive shaft. The sinusoidal velocity at the orifice was calcu- 
lated rather than measured by assuming that the sinusoidal volume 
velocity generated by the piston was equal to the sinusoidal volume 
velocity through the orifice. This relationship was assumed to be 
justified because of the stiffness of the cylinder walls and the 
relatively high incompressibility of water. Most of the measurements 
were conducted at a frequency of 22 Hz. In addition to the sinusoidal 
velocity source, a steady-state velocity was introduced by means of 
a needle valve connected to a high pressure source. 

The results of Thurston's et. al orifice impedance measurements 
are shown in Figure 3. The shapes of both the orifice resistance 
and reactance curves are very similar to that of the Ingard and Ising 's 
data (see Fig. 2) with the exception of the slight increase of reac- 
tance with orifice velocity at the higher orifice values. The agree- 
ment between the shapes of the resistance and reactance curves shown 
in Figures 2 and 3 is very impressive and lends added plausibility 
to Thurston's et. al assumption that for water the orifice velocity 
may be deduced by using the law of conservation of the sinusoidal 
volume velocity rather than by direct measurement. 

2. MODEL OF THE ACOUSTICAL BEHAVIOR OF SMALL ORIFICES 

A new fluid mechanical model of the acoustic behavior of 
small orifices is described below. The model is new because it differs 
in three fundamental ways from the earlier models of Sirignano and 
Zinn. The new assumptions are: 

(1) The sound field incident to the orifice is assumed to be 
spherical rather than one dimensional as assumed by Sirignano and 
Zinn. This assumption provides the mechanism to connect unambiguously 
the relationship between the driving pressure incident to the orifice 
and the magnitude and relative phase of the orifice velocity. The 
models of Sirignano and Zinn proved deficient in this respect. For 
this assumption to be valid the orifice diameter must be very small 
relative to the incident sound wave length. 

(2) Large changes in the magnitude and phase of the acoustic 
quantities are assumed to occur in the immediate neighborhood of 
the orifice. The hot-wire measurements conducted by Ingard and 
Ising of the acoustic velocity in the immediate neighborhood of the 
orifice show that it decreased significantly from its value at the 
orifice inlet within a distance of about two to three orifice 
diameters. This suggests immediately from the theory of ideal point 
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acoustic sources that the flow near the orifice should behave as an 
unsteady incompressible flow field (providing of course, that the 
velocity at the orifice outlet is small relative to the local speed 
of soundi 

(3) The relationship between the sound pressure incident to 
the orifice and the resulting orifice outlet velocity is scaled 
according to whether the orifice impedance is constant (lf<<^ or 
varies linearly with orifice velocity (R^>X) . The orifice impedance 
measurements conducted by Ingard and Ising and Thurston, et. al 
provide this relationship directly. 

2 . 1 Approach 


The analysis starts with the equations describing the conserva- 
tion of mass and momentum written in spherical coordinates where 
only a radial inward flow u is assumed. For orifices small compared 
to the incident sound wavelength it is logical to assume that the 
flow approaches the orifice primarily in a spherical manner. The 
origin of the coordinate system is assumed to be located somewhere 
in the orifice interior as sketched in Figure 4. Assuming spherical 
symmetry, the flow field incident to the orifice will be assumed to 
be independent of the azimuthal angle (|) (defined in Figure 4a). The 
flow field contains a uniform steady-state part and an oscillating 
acoustic part. A key element of the proposed model is the use of 
the experimental data of Thurston et. al and Ingard and Ising to 
normalize the equations describing the conservation of mass and radial 
momentum of the oscillating flow field. The experimental data shows 
that two distinct regimes exist; the regimes are defined by the 
relationships that exist between the amplitudes of the incident driving 
acoustic pressure F* and the resulting amplitude of the acoustic 
velocity V*in the orifice. For sufficiently low values of P‘, the 
data showed that 


C3) 


and for sufficiently high values of P* 


p«-~(v*y 


where (_)* denotes that the term within the brackets is dimensional. 
The regime whose P*,V* relationship is characterized by Eqn. (3) is 
often called the linear regime while the regime characterized by 
Eqn. (4) is called the nonlinear regime. 
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It is clear from dimensional analysis that the proportionality 
term in Eqn. (3) must have dimensions of density times velocity 
while that of Eqn. (4) must be density. The only density term 
suggested by the physics of the flow is p*, the mean fluid density. 
Thus, for orifices exposed to intense sound levels, the nonlinear 
case described above, the relationship characterizing P* and V* is 


p*. « (yoy C4a) 


At low sound pressure levels both Ingard and Ising's and 
Thurston’s et. al. measurements show that the orifice impedance is 
dominated by the reactance X (see Figures 2 and 3) , where X may be 
written (see Eqn. 2b) 




Since I X* I >> [R* | (and thus P*/V* = Z*=X* for most practical appli- 
cations), where R* is the orifice acoustical resistance, P* is 
related to V* as follows 


(D’‘+L*)V* 


(3a) 


where D* and L* represent the orifice diameter and thickness 
respectively. Equation (3a) can also be deduced from a dimensional 
argument. If viscosity plays a negligible role in affecting 
orifice reactance (suggested by the success of Rayleigh’s slug mass 
model in predicting orifice reactance) then the only other availa- 
ble combination is o)*(D*+L*). 

Equation (3a) suggests that at low sound pressure levels 
where P* = Z*V*==X*V* that 

X*= p’^o)’^(D^-»-L^) (5) 


Figure 5 shows that orifice reactance does indeed behave according 
to Eqn. (5); the experimental data of Ingard and Ising and Thurston 
et. al. collapses into a single correlation curve. This agreement 
is remarkable when one considers the vast differences between 
these experiments. Ingard and Ising ’s measurements were conducted 
with an orifice diameter of 0.7 cm. in air exposed to sound fre- 
quency of 150 Hz while Thurston’s et. al . measurements were conducted 


6 



with an orifice diameter of 0.3 cm. in water exposed to sound 
frequency of 22 Hz. Figure 5 shows that when V* = o)*(D*+L*), 
orifice reactance and resistance are (roughly) equal. Further, 
when (D*+L*) , then R*<<X*, the so-called linear regime, 

and when V*>> 0 )* (D*+L*) , then R*>>X*, the nonlinear regime. We 
will show below that the ratio V*/o3* (D*+L*) is an important para- 
meter in the development of the analytical model, and provides a 
means of separating the two regimes. 

The successful correlation shown in Figure 5 of otherwise 
very dissimilar data suggest that P*, V*, D*+L*, are the 

appropriate quantities that characterize the changes to the flow 
field near the orifice due to the sound. Recalling that the sound 
field incident to the orifice is assumed to be spherical, the 
solution to the (spherical) equations of motion governing the 
conservation of mass and radial momentum can be vastly simplified 
by proper scaling of the various terms and retaining only those 
of importance. The idea here is to try to anticipate the order 
of magnitude of the changes of the various terms in order to 
properly normalize them. If done correctly (and here the above 
correlation serves as a guide) , then all of the dimensionless 
terms are of order unity and therefore may be rank-ordered in 
terms of their relative importance by the relative magnitude of 
their coefficients. 


To start with, we will assume that changes in acoustic den- 
sity are adabatically related to changes in acoustic pressure. 
Thus we write The other quantities will be non- 

dimensionalized lett'ing 


P* characterize the acoustic 

P* M It II 

Y* II It II 

(D=^ + L*) ** ” 

(lii*) »» " 


pressure change near the orifice 

density ” ” ” 

velocity ’’ ** " 

length scale ” 

time scale ” " 


By defining the characteristic length as (D*+L*) rather than D*, 
the effects of finite orifice thickness will be included in the 
model. In the derivation that follows, we assume that L*<<D*. 

Now introduce the nondimensional variables r, t, u, p, p. 


The last inequality on the R.H.S. of Eqn. (6) follows from the 
adiabatic relationship p* = (c*)^p*. Substituting Eqn. (6) into 
the conservation equations yields the following, the details of 
which are in the Appendix, 
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The importance of the various terms in Eqns. (7), (8), 
and C9) are determined solely by the magnitude of their coeffi- 
cients (recall that the non-dimensional terms have been normalized 
to be of order unity). To rank order these terms, the model is 
divided into the linear regime corresponding to low values of 
sound amplitude where P* = p*(o*(D*+L*) V* and the nonlinear regime cor- 
responding to high amplitudes where P* = p*V*^ . From Figure 5, 
we see that at low values of P*, R*<<X* and V*<<m* (D*+L*) . Con- 
versely, at high values of P*, R*>>X* and V*>>oj* (D*+L*) . It 
follows from our initial assumption that the orifice diameter is 
very much less than the incident acoustic wavelength (D*+L*<<X*) , 
that therefore oj* (D*+L*) >(<c* . Further, since the amplitude of 
the acoustic pressure will always be very much less than the 
ambient static pressure, it also follows that V*<<c*. 

Based on the above information, we can now rank order the 
relative importance of the various terms in Eqns. (7), (8), and 
(9) for the two regimes. 

To simplify the equations, the nondimensional parameters 
e, M, and 6 are introduced where 
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( 10 ) 


Note that M<<e or 3. The simplified equations corresponding to 
the linear regime are 


(r»u)+ (pur»).o 




diA 


-4- 


6u 






VISCOUS 



terms 


( 11 ) 


1 

2-v* 

6 




( 12 ) 


An examination o£ the various terms in the continuity 
equation (Eqn. 10) shows that (1) the second term, which repre- 
sents the divergence o£ the volume velocity, is of order unity 
and is by far the largest term, (2) the first term, which is a 
measure of the compressibility of the fluid, is very much smaller 
than the second term (recall that M<<e) and (3) the third term 
which represents nonlinear effects is the smallest. What is of 
interest here is that even for the (linear) case of small ampli- 
tude sound approaching an orifice, the flow behaves predominately 
as if it were incompressible. 

It may be of value to offer the following more physical 
interpretation of the incompressible behavior of the sound field 
in the orifice neighborhood. Ingard and Ising’s hot-wire mea- 
surements have shown that the amplitude of the acoustic velocity 
incident to the orifice increases dramatically from a very small 
value at a distance of about three orifice diameters to a relative- 
ly high value at the orifice outlet. Since this increase occurs 
over a distance very much smaller than the sound wavelength, it 
must be a hydrodynamic change rather than an acoustic change. 

To support this, assume that the length scale (D*+L*) characterizing 
the distance over which the incident (acoustic) velocity change 
occurs is equal to the sound wavelength (X*) divided by 2 tt (recall 
that time was scaled with (ca*) rather than the orifice diameter) 
Setting (D*+L*) = X*/2ir , then M = e and from Eqn. (10) 
the first term is equal to unity. The third term which represents 
nonlinear propagation effects is of order e. Thus, the resulting 
lowest order equations reduce to the classical spherical wave 
equations as one would expect. We have shown that compressibi- 
lity became important only when o)*(D*+L*) c*. This means that 

for practical (wall treatment) sized orifices (e.g. 0.15 cm. 
diameter), frequencies of the order of 30,000 Hz or higher are 
required before compressibility effects must be included. 

For the nonlinear regime, the simplified equations are 
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( 13 ) 







(|-M^P)|| 
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(14) 


(15) 


Comparing Eqn. C13) with Eqn. (10), we note that they are 
quite similar. For both regimes the flow incident to the orifice 
is essentially incompressible. Yet, in contrast to the low pres- 
sure case, the compressibility term of Eqn. (13) is smaller (by 
the fraction g) than the nonlinear term. Thus, the normalization 
for the two regimes reveals an ordering among the various terms 
that is consistent and, therefore, lends credence to the validity 
of the model. 

2.1.1 Limitations of Approach 

There are several restrictions to the proposed model which 
impose some limitations on its application. First, it is obvious 
that the flow does not approach the orifice perfectly spherically, 
but instead probably has a streaklike pattern (i.e., an instan- 
taneous flow pattern) somewhat like the dashed curve shown in 
Figure 4(b). We recognize that the instantaneous flow incident 
to the orifice is only approximately spherical- -deviations from 
a truly spherical flow are required in order to permit the flow 
to enter the orifice in an axial manner. The spherical flow 
field is singular at the virtual origin r = o. To avoid this sin- 
gularity, the spherical flow field is truncated at a hemis- 
pherical surface of radius r* = Lg defined such that at this sur- 
face the radial acoustic velocity is equal to the maximum value 
of the actual acoustic velocity through the orifice. 

Second, for all but the lowest sound pressure levels, the 
solution to the equations of motion are valid only when the sound 
is approaching the orifice- -it is not valid during the other 
half of the cycle when the sound is moving away from the orifice. 
This limitation is imposed because it has been observed by Labate 
and Ingard, and Ingard and Ising and others that even for moder- 
ately intense sound fields, the acoustic near field velocity 
separates at the orifice forming a jet-like outflow. Thus the 
flow is spherical only when it approaches the orifices except for 
the case of vanishingly small sound fields where separation does 
not occur. This case, however, has no practical interest because 
of the extremely low sound pressures involved and will not be 
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pursued. The third and final limitation is that the orifice 
characteristic acoustic velocity oi* (D*+L*) <<c* to maintain the 
incompressibility of the near field sound. 

2.1.2 Boundary Conditions 

There are three boundary conditions that must be satis- 
fied. One is that the (radial) velocity li* must vanish along 
the walls 0 = ir/2 (see Figure 4a) due to the viscous no-slip 
condition. The second is that the acoustic pressure p* must 
merge smoothly (asymptotically) into the harmonically oscillating 
driving pressure incident to the orifice. The third and final 
boundary condition is that at the hemispherical surface r* “ ^o’ 
where the spherical inflow is truncated, the acoustic pressure 
p* must be equal to the orifice back pressure. The orifice back 
pressure is always constant and equal to its ambient value (this 
is strictly true only when the orifice exit velocity is subsonic) . 
The acoustic pressure vanishes at r* =L* , the equivalent orifice 

exit. To be precise, p* = o at r* = L*°only when the flow Rey- 
nolds number L*)V'^** is sufficiently high that the flow 

separates at the orifice forming a jet. For extremely low Rey- 
nolds numbers, the whole flow field is controlled by viscous forces. 
In this case it is reasonable to assume that the flow field is 
spherical and not separated on both sides of the orifice. 


Expressed mathematically, the boundary conditions are 

written 


U 



(16) 


T ^(r,0,t)= cos(t) 


LIMI 

r — ^ oo 


r = 


L 


(D’^+Ll") 


, ©,t 


= o 


(17) 

(18) 


The model further assumes that 10 (D*+L*) <<X* so that changes 
of sound amplitudes and phase are negligible over this distance. 

2 . 2 Linear Regime 

The solution to Eqns. (10), (11), and (12) subject to the 
boundary conditions defined by Eqns. (16), (17), and (18) describe 
the behavior of orifices exposed to weak to moderate sound fields. 
Retaining terms to order e (since NJ << e, the effects of compres- 
sibility would only weakly affect orifice impedance and are thus 
ignored) yields 
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C19) 


^(r^u) = 0 

^|> _ _2 J_ ^ _ Q 

^0 Re r 3© 


( 20 ) 


( 21 ) 


where Re = to* (D*+L*)^/v* is the orifice characteristic Reynolds 
number. Equation (19), representing the conservation of mass, 
shows clearly that the flow oscillates through the orifice in an 
imcompressible (and unsteady manner. 

Much can be learned about the behavior of orifices by 
carefully examining Eqn. (20). To begin with, assume both 
e = 0 and v = 0 and further that both u and p are harmonically 
driven so that u(r,-t) = d (>") and p^r,-t)= e'* p(*~) 

Under these circumstances it is clear that u is ninety degrees 
out of phase with p and that the orifice impedance is reactive 
only (which it should be if the fluid is frictionless) . Now 
allow the fluid to be real (i. e., viscous). For most practical 
applications, the orifice characteristic Reynolds number is such 
that a laminar boundary layer forms along the walls containing 
the orifice as sketched in Figure 4(b) (assuming typical values of 
f=1000Hz, (D*+L*) =0 . 15cm, then Re~10^ which is assumed to be in the 
laminar Reynolds number range) . Thus the resistive losses due 
to viscosity should be small and hence orifice resistance should 
be very much less than reactance. Now consider finite values of 
e. From Eqn. (20), it is clear that the loss of the fluid kinetic 
energy e pu^/2 is in phase with the pressure (i.e., resistive) 
and increases with increasing u. In summary, at very low values 
of sound pressure, orifice reactance is very much larger than 
resistance; with increasing sound levels, the resistance should 
increase because of jet kinetic energy-type resistive losses. 

This is in complete agreement with the measurements of Ingard 
and Ising and Thurston et. al. (see Figure 5). 

The formal solution to Eqns. (10), (11), and (12) consist 
of expanding u and p in powers of e 


u(r,e,-t) = u,(r,Q,t) + £ u, (r, ©,t ) 4- • • • (22) 
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+ 6 f , (r,0,t)+ ••• 


(23) 


Substituting Eqns. (21a) and (21b) into Eqns. (10), (11) and 
(12) yields to order e° and e 


and 




+ — fj- 4 - (r 

at ar cO^(d^+LT)*- L r*- ^ V 

q ! 4- (sino -1^)1 

sin© a© V 3© ! 


% 

bv! 

= o 


2u^ 

r*- 


3^;; 2V* _!_ 

be uJ*(\>*+lf)^ ^ 3© 


(24) 


(25) 

(26) 


A 

3r 



= o 


au, 

3i: 


+ 




t;* 

u)*(D*+L*)*' 


V/ISCOUS TERMS 


= o 


_i_ 3^ 

a© co^CD*+L?)*' r a© 


(27) 

(28) 

(29) 


The solution to the above equations must satisfy the bound- 
ary conditions specified by Eqns. (16), (17), and (18). The usual 
procedure in seeking perturbation solutions is to force the lowest- 
ordered terms to satisfy the boundary conditions with the higher- 
ordered terms set equal to zero. Thus the appropriate non-dimen- 
sionalized boundary conditions are 
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and 


and 


u.(r,© 

LIMIT ■R( r,©,-t-) = cos(-t) 

r ->oo ^ / 


l: 


lT) 


) o, t 


= o 








r = 


L. 


(D* + U*') 


jOj't 


= 0 


(30) 

(31) 

(32) 


(33.) 

(34) 

(35) 


2.2.1 Lowest Order Solution 


The solution to Eqns. (24), (25), and (26) subject to the 
boundary conditions given by Eqns. (16) and (17) represent a first 
approximation to the behavior of the acoustic flow near the orifice. 
For convenience, the solution will be divided into two parts. In 
part 1, the effects of viscosity will be ignored. This means that 
the no-slip boundary condition given by Eqn. (16) will be ignored. 

The solution to part 2 will include the effects of viscosity. 

The solution to Eqn. (24) applies equally to both the inviscid 
(part 1) and viscous (part 2) solutions and may be immediately 
integrated to yield 


F^(©,-t3Re) (36) 

where A is an arbitrary constant and F is an unspecified function. 
The negative sign is included to denote that the velocity is 
directed radially inwards. 
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Part I - Inviscid Solution 

The equations describing the inviscid behavior o£ the fluid 
near the orifice follow from Eqns. (25) and (26) by setting v*=0 
to yield 







(37) 

(38) 


Equation (38) shows that Pq is independent o£ the azimuthal angle 
6. Since p^ is independent of 0, it also follows from Eqn. (37) 
that Uq is independent of 0 and hence from equation (36) that Fq 
is also independent of 9. Thus, equation (36) becomes 


Substituting Eqn. (39) into Eqn. (37) yields 





r cl-t 


= o 


which integrates to 


Cnt) 


A dF. 

r d-t 


= ■fo('t) 


(40) 


where £ (t) is an arbitrary function of time. From the boundary 
condition given by Eqn. (17), f (t) - cost and Eqn. (40) may be 

written 

9.(r,-t)= cos(t)-A (41) 


According to our model, the radial inflow accelerates as it 
approaches the virtual source, hence it must at some radial distance 
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reach a magnitude at which it is equal to the magnitude of the 
maximum orifice outlet velocity. It follows since the outlet 
flow is subsonic that the outlet pressure must be equal to the 
ambient background pressure. Thus the acoustic perturbation at 
the orifice outlet must vanish which means that 


r = 


CD^+U^) 


A = L.V, 


(D*+Lf^) 


( 42 ) 


Since the acoustic velocity has been normalized to unity at the 
orifice outlet, we see from Eqns. (39) and (42) that A = r = 1. 
Thus, the velocity and the pressure written in complex notation 
for convenience are 

u„(r,-t)= -linM = (43) 

“V ) ; pi 

and 

W = ('- -f) 


Equations (43) and (44) are our desired results. Following Ingard 
and Ising, we define the orifice impedance as the ratio, of the 
driving pressure far (in terms of the characteristic length 
(D*+L*) from the orifice to the velocity at the orifice outlet. 
Dividing the pressure Po = far upstream of the orifice by the 

orifice outlet velocity = ie yields for the orifice 

impedance 


e(r=l,-t) 


In dimensional terms, the impedance is 




(46) 


There is no resistance, only a negative reactance, because 
the fluid has been assumed to be inviscid (in the linear sound 
regime the only mechanism to dissipate acoustic energy is through 
the fluid viscosity) . The negative sign occurs because the inward 
travelling sound is negative in a spherical coordinate system. 

Both the orifice resistance and reactance for the linear and 
nonlinear regimes will be negative for this reason. The orifice 
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reactance given by Eqn. (45) agrees exceptionally well with the 
experimental measurements taken by Ingard and Ising and by Thur- 
ston, et. al. which is shown in Figure 5. The reactance derived 
above is virtually identical with the reactance derived by the 
slug-mass model (Eqn. 2b) , thus showing the equivalence of the 
two models at low sound pressure levels. 

Part 2 - Viscous Solution 


At low sound pressure levels, Ingard and Ising 's measure- 
ments showed that the (linear) resistance of an isolated orifice 
may be approximated by the empirical expression repeated here as, 

R*s \J 8 jU*e*co* (| + >-*/D’)i X*S(oW(u*+0.85D') ( 


Ingard^ and Crandall’s^^ work led to the derivation of Eqn. (2). 

They solved for the relationship between the pressure gradient, 
unsteady velocity, and viscosity in an infinitely long cylindrical 
tube containing an unsteady, viscous, fully-developed one-dimen- 
sional flow. Ingard^ later corrected this to include finite ori- 
fice thickness by including as an end correction, the viscous 
contributions of the side walls containing the orifice. This lead 
to the development of Eqn. (2). A major deficiency of this approach 
is that for most practical orifices of interest, the orifice thick- 
ness is small compared to the orifice diameter rather than very 
large as the derivation requires. In this sense, the end correc- 
tion is not a correction at all but represents instead the major 
part of the losses (i.e., the losses associated with the viscous 
flow along the side walls containing the orifice) with the losses 
inside the orifice being small. 

For orifices where L*<<D*, the use of an end correction to 
account for viscous losses as given by Eqn. (2a) is unnecessary. 

It follows logically from the proposed model that the viscous 
losses arise from the boundary layer established by the radial 
flow moving inwards along the side walls containing the orifice. 

The importance of viscous effects are characterized by the flow 
field Reynolds number. If the Reynolds number is large, then 
viscous effects are important usually within a small region called 
the boundary- layer region. Conversely, if the Reynolds number is 
small , viscous effects are important throughout the whole region 
of interest. Assuming values representative of the orifices used 
as acoustic liners in present day jet aircraft, that is 

10* HZ., 1^^-0.15‘cm^ec, then the Reynolds number Re = 10* 
which is sufficiently large to model viscous effects using 
boundary -layer theory. 
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Boundary -layer theory assumes (1) that viscous effects are 
important only within a small region, in the present application, 
adjacent to the side walls containing the orifice, (2) the pressure 
gradient normal to the wall is negligible, hence the pressure is 
constant through the boundary layer and its value is derived from 
inviscid theory and (3) that the velocity at the wall surface is 
zero and approaches asymptotically at the boundary layer edge its 
inviscid value. The appropriate boundary- layer equations to be 
solved are 






C47) 



Ao - _L 

— — i 

auA_2uo 1 


sin© 

at 

ap R.e 

ar ' 

’*’r*sine 

a© ^ 

be) 


=0 (48) 


^0 


= o 


(49) 


Equation (48) can be simplified further by introducing the 
transformation X €> and observing that since the boundary 

layer is confined to the immediate neighborhood of the wall where 
B = V</'Z Figure 4b), then ^ is very small and we may appro- 
ximate 5ir\0 = cos^5) “ 1 . Introducing the angle ^ and the above 

approximation into Eqn. (48) yields 


at ar 



J_ A (r^ - Ai + _L 

arV 


= O 


(50) 


Eqn. (47) may be integrated and inserted into Eqn. (50) which, 
after some algebraic manipulation, yields 


u,(r, Re) = 


ie'^ 


(51) 


where we have used complex notation and note explicitly that the 
function Fq depends both upon the angle ? and the Reynolds number 
R . From boundary-layer theory and from the no-slip boundary 
condition, Eqn. (16), Fq(?. Re) must satisfy that 

Fo(^= 0 )Re) = 0 and ^e) = 1 ( 52 ) 

t;yo 
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The boundary condition on the function F has been selected 
so that the viscous velocity u approaches asymptotically (with 
increasing Reynolds number) its inviscid value defined by Eqn. (43). 
Using Eqn. (44) for the (inviscid) pressure and substituting u^ 
and Po into Eqn. (50) yields after some algebra 




Re d ; 




(53) 


To recover the boundary- layer character of Eqn. (53), we introduce 
the boundary- layer coordinate n defined by 

y r ^ (54) 


and set F (^,Re) = * Physically, this means that the 

velocity is ?nvariant to the transformation defined by Eqn. (47). 
Substituting Eqn. (54) into Eqn. (53) yields 

I - =0 CSS) 


The solution to Eqn. (55) that satisfies the boundary conditions 
given by Eqn. (52) is 

-(l+i) s/T y 

h(y) = I - e ^ ( 56 ) 


Substituting Eqn. (56) into Eqn. (51) and solving for the real 
part yields for the velocity Uq 




■(n7.-t)= -pr 




( 57 ) 


The instantaneous velocity distribution (Eqn. 57) is uniform 
everywhere near the orifice except near the walls 6 = tt/2 where 
the velocity decreases from its uniform value to zero through a 
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thin (boundary layer) region of thickness Re^l^ . The retarding 
action of the fluid viscosity acts, for a given driving excitation 
pressure, to decrease the magnitude of the acoustic velocity pumped 
through the orifice by an amount related to the displacement 
thickness 6* defined as 


U, 






Substituting Eqn. (57) into Eqn. (58) yields 




(58) 


(59) 


Since Uq 6* represents the amount of mass "lost" because of the 
effects of viscosity, the effective velocity approaching the 
orifice is 


U, 




_ !_^ _ COsW 


(60) 


At r = 1 (the effective orifice outlet), the velocity is, written 
in complex notation, 


U 


0 ,+) 


= e'^ 


iZRe I {2^ 


(61) 


Dividing the pressure p = e^^ incident to the orifice by the 
velocity Uq given by Eqn. (61) yields for the orifice impedance 


(-z) s 




4- 


Re 


+ 


1 + 


\j2 Re 


( 62 ) 


Resolving Z into its real (I'^sis tance) and imaginary (reactance) 
parts yields 
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and 


( 63 ) 



(64) 


Written in dimensional terms, the impedance is 



4- 




(D^+L*)J 




(65) 


Equation (65) is o£ interest because it shows the behavior o£ ori£ice 
resistance as namely (-R)^y*/(D*+L*) . Thus, we see that the 

e££ects o£ viscosity are to generate a resistance given by Eqn, (63) 
and to perturb the inviscid reactance (Eqn. 45) by the small (con- 
stant) amount given by Eqn. (64). 

The ori£ice (linear) resistance was derived using the concept 
o£ a boundary- layer displacement thickness which may be inter- 
preted physically as a loss o£ mass £low along the walls contain- 
ing the ori£ice due to the retarding e££ects o£ viscosity. The 
derivation demonstrates clearly that the proposed spherical model, 
based on £luid mechanical concepts, is the proper way to understand 
the acoustic behavior o£ ori£ices. In contrast, Ingard and Ising 
interpret their (linear) ori£ice resistance data (£or their ori£ice 
L*<<D*--see Figure 1) in terms o£ a *Viscous end correction”. But 
the £undamental derivation o£ this ”end correction” by Rayleigh and 
others assumed that to £irst order the ori£ice thickness was 
in£initely long (i.e., L*>>D*) . Their derivation assumed that the 
oscillatory motion o£ the £luid contained within the ori£ice 
thickness was equivalent to a solid body oscillation. To £orce 
agreement between theory and measurements an ”end correction” was 
added. To apply this model to predict the resistance o£ their 
ori£ices where L’^<D*is clearly unsatis£actory although the empi- 
rical £ormula (Eqn. 2) is in agreement with experimental £indings. 
This approach, however, breaks down when it is extended to predict 
ori£ice nonlinear resistance. The solid body oscillation approach 
is clearly incompatible with the observed nonlinear jetting or 
•’breaking away” o£ the £luid £rom the ori£ice. In contrast, the 
model presented herein does not have these de£iciencies , but in- 
stead includes, in a straight £orward manner, the important e££ects 
o£ viscosity and jetting on the ori£ice resistance and reactance . 

The predicted linear ori£ice resistance (Eqn. 63) di££ers 
£rom the empirical prediction (Eqn. 2a) by the factor (2). Apply- 
ing both equations to the orifice viscous resistance data measured by 
Ingard and Ising (summarized in Figure 2--here L*<<D^ , it would 
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appear, at first glance, that Eqn. (2a) was more accurate. We 
shall show that this may not be true because nonlinear effects 
will be shown below to be important even at the low orifice velo- 
cities . 

2.2.2 First-Order Solution 


The differential equations v (27) , (28), and (29) represent 
nonlinear corrections (of order ) to the lower ordered 

linear solution. To simplify the analysis, only the inviscid 
correction will be sought. Thus,T?*^= o and the equations describing 
the first order corrections become 



II 

o 

(66) 


4- ^ = O 

or or 



(67) 

Ai = 

a© 

0 — ?, =P, (r, t) 

(68) 


Equation (67), in its present form, is misleading. By 
rewriting it in the form 



it is clear that the term 1/2 uq^ is in phase with the pressure p^. 
Since it is always positive, it must be replaced with 1/2 Uq|uo 1 
if it is to be in phase with the harmonically fluctuating pressure 
This restriction is based mostly on Ingard and Ising's hot- 
wire measurements which showed that u^ fluctuated harmonically (at 
least along the orifice center line) . ^An alternate explanation is 
that the model is valid only during the half-cycle during which 
the acoustic velocity is approaching the orifice; to account for 
the other half, the coordinate system must be mathematically 
"switched*’ to the other side of the orifice. While the linear 
terms account for this automatically, the non-linear term 1/2 Uq^ 
does not. It follows from this discussion that the correct momen- 
tum equation is 


^ + A 

at ar 




(67a) 


Integrating Eqn. (66) yields for the first-order velocity 




( 69 ) 
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where Fi (t) is an arbitrary function of time. 

Substituting Eqn. (69) into Eqn. (67a) and integrating with respect 
to r yields 

where fj (t) is also an arbitrary function of time and u^ is given 
by Eqn. (43). From the boundary contition specified by°Eqn. (33) 

Limit ^(pt) = 0 ->f, (-t) = 0 

r ->oo 

Thus, Eqn. (70) reduces to 

-L + 

For application to nonlinear flow, the model assumes that at the 
radial location r = 1, the flow + eui is accelerated to its 
maximum value and the pressure p^ + epi OP^'oviding the maximum 

velocity u + eui is subsonic) is equal to the ambient background 

pressure. °Thus, the acoustic perturbation px must vanish when the 
velocity is a maximum and we write at r = 1 


(71) 


(72) 


f^(r=l,t)=o 


Substituting Eqn. (73) into (72) yields at r = 1 


a-b 


- X I X lsin(t)lsin(t) 


(73) 


(74) 


The integration of Eqn. (74) ds split into two parts, one corres- 
ponding to t>o and the other to t<o. Thus, Eqn. (74) integrates 
to 



sin(t) cos(-fc) 

J ) 


t< o 


sin(t)cos(-fc) 


-t > O 


( 75 ) 
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where the constant o£ integration has been ignored. 


Since we are interested primarily in predicting orifice 
impedance, only the fundamental harmonic component of Eqn. (75) 
is of interest. Using Fourier analysis, and noting that Fi (t) 
is an even function, it is straight-forward to show that 


a, = :| J^F;(-t)cos(^:)cHr = '^Tr 

o 

and hence that 


( 76 ) 


F;(t)sjL tOs(+) (77, 

Combining Eqns. (69) and (77) yields for the nonlinear velocity ui , 

(78) 

The orifice velocity, including both viscous and nonlinear contri- 
butions, follows (written in complex notation) by combining Eqns. 

(60) and (70) to yield. 


U 



pa 




f 


37T 


V2Re 


■) 


(79) 


Again, defining orifice impedance as the ratio of the 
driving acoustic pressure to the orifice velocity, it follows that 






The real and imaginary parts of Z are 



(80) 


(81) 

( 82 ) 


Both Ingard and Ising and Thurston et. al presented their 
orifice impedance data in terms of the ratio of the incident 
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driving acoustic pressure to the orifice inlet velocity. The 
definition of impedance used in this model is based o/n the maximum 
orifice velocity rather than the inlet velocity. Before the orifice 
resistance and reactance defined by Eqns. ( 81 ) and ( 82 ) can be 
compared with the experimental data of Ingard and Ising and Thur- 
ston, et. al, these differences must be resolved. Significant 
differences occur between the orifice inlet and maximum (vena 
contracta) velocities only for the part of the orificq velocity 
that is in-phase (resistive) with the incident driving pressures. 
Orifice flows that are in-phase with pressure can be characterized 
as quasi-steady and the concept of a discharge coefficient Cj, 
(suitably time -averaged) can be introduced to connect the orifice 
inlet velocity V. to its maximum value V (at the vena contracta) 
where ^ 


Vi 


= 


(83) 


In steady-state flows through orifices, the discharge 
coefficient is used because it conveniently determines the mass 
flow through the orifices for a given driving pressure. In most 
applications, the characteristic flow Reynolds number is suffi- 
ciently high (above about 1000) that the discharge coefficient is 
constant. The solution of the nonlinear acoustic velocity ui 
shows that (1) it is in phase with the driving pressure p (see 
Eqn. 78), (2) it is quadrat ically related to the u (see Eqns. 44 and 67) 
and (3) it is independent of viscosity (this corr§sponds to flows 
of infinite Reynolds numbers) . It is clear that Ui is quasi- 
steady and inviscid and falls within the category of a Bernouilli- 
type of flow. Thus, to connect the model results with experimental 
data, Eqn. (83) is used to replace V by V. in the definition of the 
perturbation parameter e rewritten as 


£ = 


V: 




(84) 


Substituting Eqns. (83) and (84) into Eqns. (81) and (82) yields the 
final form of the resistance and reactance. 


( ViRe R«.) 3n c-o L"* 


V;^ 


(i 




Re ' Re') 


(85) 


and 



r 

1 ® 1 

1 


y 1 

3TT Cp 1 




( 86 ) 
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Transforming Eqns. (85) and (86) to dimensional terms 



viscous loss Je-t loss 

and 


( 87 ) 


V '' -2/ 

OScillatiri^ slu^ mass 
model recovered 


6 Vi* 

37V oo^(^0*+L^) 


(88) 


The negative signs occur in Eqns. (87) and (88) because 
the sound field is incoming. The orifice resistance described by 
Eqn. (87) consists primarily of the sum of two terms, one repre- 
senting a nonlinear Bernouilli-type of jetting effect and the other 
a viscous dissipation loss effect. The orifice reactance, described 
by Eqn. (88), is also the sum of two terms, one representing an 
oscillating mass flow term and the other a higher-ordered combined 
viscous and nonlinear (Bernouilli) term. Equation (88) shows the 
nature of the cross-coupling between the nonlinear Bernouilli 
effect and orifice reactance. The effect is to weakly decrease 
orifice reactance. The effect is weak because the derivation re- 
quires that Vj^*<<o) (D*+L*) which restricts the magnitude of the 
jetting effect. The existence of a weak reduction in orifice reac- 
tance has been observed experimentally over a wide range of test 
conditions including the experiments by Ingard and Ising with air 
(Fig. 2) and Thurston, et . al . with water (Fig. 3). A direct com- 
parison between the reactance predicted by Eqn. (88) and the ex- 
periments by Ingard and Ising and Thurston, et. al is shown in Fig. 5. 
The comparison is excellent. 


Since Vi^co*(D^+L^)« I , it follows that |R|<<|x 1 and 
that the orifice impedance ^ -X is virtually constant. 

Although the nonlinear Bernouilli term only weakly affects the 
reactance and hence the impedance, it can dominate the viscous 
resistance term in Eqn. (85) providing that 


3TT 

^\jZRe 


« 


Vj* 

(P*+L^) 


« t 


( 89 ) 


Thus an interesting behavior or embedding exists wherein the resis- 
tance is dominated by nonlinear (jet) losses of the kind described 
by Ingard and Ising but the impedance is dominated by an almost 
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constant reactance - hence the impedance is, to first order, constant. 

To verify the usefulness of Eqn. (87), it should be compared 
with the experimental data of Ingard and Thurston. Unfortunately, 
the behavior of the discharge coefficient at these low velocities 
is not available. As an indirect means of verifying its usefulness, 
Eqn. (87) will be forced to agree with Ingard* s and Thurston's 
measurements by computing the required values of Cj. The results 
of this force fit are shown in Fig. 6 in terms of the required 
values of Cj.. The shape of the required Cp curve is similar to the 
measured Cjj aata shown in Melling's paper, which shows the effect 
of Reynolds number. The similarity between these curves despite 
their difference in diameter (D*) , thickness to diameter ratio 
(L*/D*) , and the fact that one represents a force fit to a.c. 
resistance data and the other d.c. test data is encouraging. To 
verify the usefulness of the model and hence Eqn. (87 j, detailed 
a.c. tests of the dependence of Cj) to geometry and Reynolds number 
are required. 

2.3 Nonlinear Regime 

The basic equations governing the behavior of intense 
sound fields near orifices were derived in Section 2.1 above and 
are defined by Eqns. (13), (14), and (15). In the nonlinear regime, 
the characteristic orifice Reynolds number. Re = V*(D*+L*)/v* 
is independent of the frequency in contrast to the Reynolds number 
characterizing the linear regime. Re = w*(D*+L*) ^/v*. Since the 
Reynolds number is fairly large in the nonlinear regime (typical 
jet engine linear values are V* = 50 m/sec, (D*+L*) = 1.5x10 “®m, 

V* = 0 .15x10 ~‘*m^ /sec to yield Re ^ 5000), the contribution of 
viscous terms to the orifice impedance will be ignored. 

Retaining terms to order 8 and ignoring viscous terms, 

Eqns. (13), (14), and (15) simplify to 
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II 

(90) 



(91) 
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o 

II 


(92) 


The velocity follows immediately from Eqn. (90) to be 

= ~ F(t)/r*- (93) 

where F(t) is an arbitrary function of time. Substituting Eqn. (93) 
into Eqn. (91) and integrating with respect to radial distance (r) 
yields 
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where g (t) is yet another arbitrary function of time. Since vis- 
cosity is ignored, the no-slip boundary condition specified by 
Eqn. (16) does not have to be satisfied. The pressure, however, 
must satisfy the boundary contitions specified by Eqns. (17) and 
(18) . Applying the boundary condition specified by Eqn. (17) to 
Eqn. (94), it is clear that g(t) = cos (t) . Applying the boundary 
condition specified by Eqn. (18), p(r = 1, t) = o, yields the fol- 
lowing equation at r = 1 


The solution to Eqn. (95) describes the time behavior of 
the vena contracta velocity (i.e., the maximum velocity of the jet 
exiting the orifice) . Equation (95) is a nonlinear differential 
equation of the Riccati type. Since, to our knowledge, there are 
no known analytical solutions to Eqn. (95) , an approximate solution 
is sought sufficiently accurate to elucidate its physical interpre- 
tation. Since 3<<1, the contribution of the unsteady term is relatively 
unimportant; the orifice behaves in an essentially quasi-steady manner. 

The usual procedure to solve Eqn. (95) in an approximate way 
is to expand the function F in powers of 3 and to solve, say, for 
the first two functions F^, Fi defined below as 

F(-t,p) = F.(t) + pF, (+)+••• C96) 


This kind of approach, however, breaks down and the solutions to 
Fp and Fi become singular for values of t = (2n - l)ir/2, n = 0, 1, 
2.... To illustrate this singular behavior, the solution to Eqn. 
(96), assuming that 3=0, is 


Fo(t) = \I 2. cos(±) 


(97) 


The derivative of (F) in Eqn. (97) is of order unity every- 
where except near t =(2n-l)ir/2 where it becomes singular. Physi- 
cally, the flow oscillates through the orifice with period 2Tro)"^ 
and near t = (2n - l)ir/2 the function Fq is both small and rapidly 
decelerating (or accelerating). Thus, near t = (2n - 1 )tt/ 2 the 
assumption that «• is incorrect and both terms 

in Eqn. (95) must be retained. 

The singular behavior of F(t, 3) is analogous to the singular 
behavior of boundary- layer type flows wherein large changes to the 
flow tangential to a surface occur within a small region. A suc- 
cessful procedure to remove such singularities is to use the method 
of singular perturbation theory. The details of the derivation that 
follows are described in the Appendix. To start with, Eqn. (95) 
is transformed to the boundary layer form 
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( 98 ) 


where 

F^jp)= ^ \fz p 

To solve Eqn. (98) , the boundary- layer type function f is expanded 
as follows 


f(4,J) - •f.(4) + S*'^"'f,(4) + o(f/3) (100) 

Substituting Eqn. (100) into Eqn, (98) and collecting the coef- 
ficients of 6° and 6 yields 

- ii q- = t (101) 

-^ + ^-if, =-i 4" 

2.3.1 Solution to Order 8° 


Equation (101) is the well-known Riccati equation whose 
solution can be expressed in terms of the Airy function Ai(£i) as 


t(4) = 


(103) 


where (*) denotes the operator d/d?. A plot of the functions 
-fo(4) aiid /y4 are shown in Figure 7. It is evident that 
the function has the correct asymptotic behavior to properly 

match to the inner behavior of (the details of the matching 

of the inner and outer functions are given in the Appendix) . De- 
fining the thickne ss o f the inner region (?-,„ ) as the value of ?»« 
where -f„ (%tO = 0.99 , then from Figure 7, 5. Thus, to 

first order the solution to Eqn. (95) is 
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It is convenient to rewrite Eqn. (104) in the following form 

F(t) s \/ 2 cos ft) + [^W] 


where Af (t) is defined in Figure 8. Substituting Eqn. (105) into 
Eqn. (93) , the radial velocity approaching the orifice is approx- 
imately 

u(r,-t) = -f(^)/^z ^ ~y^z |y^2cos(-t) 4- 

The harmonic component o£ Eqn. (106) is determined approximately 
by Fourier analyzing the resulting expression yielding 





- i.S7 cos(^) - 


sin (-t) 


2.3.2 


Solution to Order 


B 


4/3 


The solution to Eqn. (102) represents the first-order cor- 
rection to fQ(5). The homogeneous solution to Eqn. (102) is 




(104) 


(105) 


(106) 


(107) 


(108) 
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where A is an arbitrary constant of integration. The particular 
solution is 



(109) 


The general solution to Eqn. (102) is 


f,(4) ® ( 4 )+ -(-,^(4) 


( 110 ) 


The only way to match the outer solution of Eqn. (110) to the inner 
behavior of F(t;3) is to set A = 0. Thus, to order , the solution 
to Eqn. (95) is 



^ si Zcos(^) , O <-t < ^ - 5 ( 2 ^ 3 ) 


( 111 ) 


A plot of the functions fi(4) and its asymptotic behavior -1/12 4 

is shown in Figure 9. In a manner similar to that used in the de- 
rivation of the approximate solution given by Eqn (104) , a solution 
to order of Eqn, (95) is 




where the function Afi [ ?(t)] is defined in Fig. 10. The harmonic 
component of Eqn. (112) is determined by Fourier analysis to yield 


F(t)|S) = - 1.57cos(t)-(l.62p+ 


(113) 


Thus the harmonic component of the radial inflow ui written in com- 
plex notation for convenience, is approximately 



pi 


(114) 
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Following Ingard and Ising, we define orifice impedance as the 
ratio of acoustic pressure far from the orifice to the orifice inlet 
velocity V;*. As described earlier in Section 2.2.2, the relationship 
between the orifice inlet velocity and the orifice maximum velocity 
V* (vena contracta velocity) is defined as 





(115) 


where is the time-averaged discharge coefficient. Equa- 
tion (li5) is valid only for the part of the velocity that is in 
phase with the driving pressure. Replacing ui(r,t) at r = 1 with 
in Eqn. (114) yields the following estimate of the orifice 
outlet velocity u. written in complex notation 




= -157ci^e'*+ i (l.62p+4.04p''/Oe’'^ 


(116) 


Using Eqn, (116), the orifice impedance is 


(-z) = 


(117) 


1.57 Cj) l(l.fe2|j+4.o4|a^)e'* 


The orifice resistance and reactance follows from Eqn. (117) 
to be approximately 


(-R)- 


0.64 




(118) 


and 


(-X) * ^ 


+ 


2.48f 




(119) 


In dimensional terms, the resistance and reactance are 




- I 


( 120 ) 
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and 
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To compare the predicted orifice behavior with the data by Ingard and 
Ising and Thurston et.al, Eqns. (120) and (121) are nondimensionalized 
by the quantity p*o)*(D*+L*) to yield the final form 


r 

\ 

1 


1 



15) 


CO* (D*+L*) 


[' Vi- J\ 


( 122 ) 


and 
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Equations (122) and (123) must be interpreted with caution 
because of two key assumptions made in their derivation. First, the 
solution is valid only when the sound field is approaching the orifice 
it is not valid during the other half of the cycle when the sound is 
moving away from the orifice. Second, the orifice reactance was 
derived only near t = tt/2 (where the flow direct ioi changes) - it was 
assumed to be zero at t = 0 . It is clear from Ingard and Ising’ s 
measurements , that the phase shift between the driving pressure and 
the orifice velocity is constant for a constant driving sound pressure 
and frequency. Thus a reasonable correction to Eqn. 123 would be to 
reduce its value to one-half under the assumption that by forcing the 
reactance to be zero at t = 0, the reactance at t = tt/2 is likely to be 
twice its true value. This is possible even with the first assumption 
because Ingard and Ising have shown that the orifice inflow is symmet- 
rical during each half-cycle. Accordingly, Eqn. (123) is reduced by 
one-half to 
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Equations (122) and (124) show that orifice resistance and reac- 
tance are coupled through C„ the discharge coefficient. To verify 
the utility of the model, tne orifice resistance and reactance 
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measured by Ingard and Ising for air (D = 7mm, f = 150Hz) and by 
Thurston, et. al for water p = 3mm, f = 22Hz) are compared with 
the predicted values determined by Eqns. (122) and (124). Choosing 
Cjj = 0.80 the predicted values compare remarkably well (particularly 
tne shape) with the measurements for both very different test condi- 
tions (for air and for water) as shown in Fig. 5. 

The model shows clearly that at high sound pressure levels 
where V. *>>o)* (D*+L*) , orifices behave in a predominately quasi-steady 
manner In agreement with the experimental findings of Zorumski and 
Parrott^®. Both orifice resistance and reactance are shown to be 
related to a time-averaged discharge coefficient. The quasi-linear 
behavior of the orifice is now clear. Over most of the cycle the 
orifice reactance term | 8 9u/ 3t | << lu3u/3x 1 except near the turning 
points where t = (2n-l)ir/2. The analysis shows that reactance is 
important only when ir/2 - (5) (2 ^'^ ) (8^'® ) <t<T:/2 or about (5) (2 ^/®) (8 ^^) /tt/ 2 
of a cycle. Physically, this means that the slight phase shift between 
the driving pressure and the orifice velocity is unimportant over 
most of cycle but becomes important when the orifice velocity is 
sufficiently low (of the order 8*^) that the slight phase shift 
dominates and the orifice flow direction reverses. 

3. SUMMARY OF RESULTS 

The results of this model extends significantly our understand- 
ing of the acoustic behavior of orifices, particularly reactance. 

The model has direct and important application to the connection 
between the optimum wall impedance required for maximum sound attenua- 
tion and the wall construction having this desired impedance. 

The results of this study show the following: 

. The acoustic flow in the immediate neighborhood of the orifice 

(i.e., the near field) can be modelled as a locally spherical flow. 
Within this near field, the flow is, to a first approximation, un- 
steady and incompressible. This is true regardless of the intensity 
of the incident sound pressure field. Thus the behavior of orifices 
are hydrodynamic rather than acoustic. 

. The behavior of orifices can be roughly divided into three 

regimes depending upon the value of the ratio V. */aj* (D* + L*) where 
V. is the magnitude of the orifice inlet velocity, oj* is the sound 
ridian frequency, and (D*+L*) is the orifice diameter and thickness 
respectively. When V. (D*+L*) <<1 , orifice reactance is much 
larger than resistance and the orifice impedance is constant, indepen- 
dent of the incident sound pressure (the so-called "linear" regime). 
When V . */oj* (D*+L*) >>1 , orifice resistance is much larger than 
reactaAce and the orifice behaves in a quasi-steady manner, (the so- 
called "nonlinear" regime). When V- */o3* (D*+L*) = 0(1), orifice 
resistance is approximately equal to reactance. 

The parameter V. (D*+L*) can be directly related to the 
incident sound field. ^ In the linear regime, V.* = P*/p*(o* (D*+L*) 
where p* is the fluid mean density and P* is tne amplitude of the 
incident sound field. Thus, in the linear regime, the ratio 
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V|/o3* (D*+L*) * P.*/p* [oj* (D*+L*) ] ^<<1 . In the nonlinear regime, 

V.* = [CjjP*/p*] where Cj, is an orifice time -averaged discharge 
coefficierLt . Thus, in the nonlinear regime, the ratio V^/a)*(D*+L) = 
(Cj.P*/p*) (D*+L*) >>1 (to be consistent with the paraAeter used 

above for the linear regime, one could use, since Cj. 0(1), 

P*/p* [o)*^ (D*+L*) ] ^>>1 to identify the nonlinear regime). The inter- 
mediate regime where orifice resistance and reactance are roughly 
equal can be identified when P*/p* [w* (D*+L*) ] ^ 1. 


. For P*/p* [o)* (D*+L*) ] ^<<1 , corresponding to a reactance domin- 

ated orifice impedance, orifice resistance consists of two terms, 
one a viscous term related to the boundary layer displacement along 
the walls containing the orifice and the other a jetting type term 
related to the flow nonlinearity. With the exception of extremely 
low sound fields, the nonlinear resistance dominates the viscous 
resistance; its behavior is expressed in terms of a time-average 
discharge coefficient. Comparison between the model and data shows 
that the discharge coefficient is a function of the ratio 
P*/p* [oj* (D*+L*) ] ^. Despite the existence of the nonlinear resis- 
tance term, the orifice impedance is virtually constant because 

I Xl»|R| . 


For P*/p [co* (D*+L*) ] ^>>1 , corresponding to a resistance dominated 
orifice impedance, the orifice behaves in a quasi-steady manner in 
agreement with the experimental findings of Zorumski and Parrott 
Both the orifice resistance and reactance are shown to be related to 
a time-averaged discharge coefficient. 


. The results of the model have been compared with the orifice 

impedance measurements by Ingard and Ising in aiv and by Thurston 
et. al in water and the agreement between prediction and measurement 
is excellent. 
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APPENDIX 


SUPPORTING ANALYSIS 


To simplify the amount of analysis presented in Section 2, many 
of the supporting or intermediate steps were omitted. For those 
interested, they are described below. 

Equations (7) , (8) (9) 

Assuming spherical symmetry, the flow field incident to the ori- 
fice is 
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(A-l) 


■ I ha 

- 2u* I ^ /'sin© _ 

r* r»*^sin© a© \ b® > ^ 


J ^ f 
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CA-3) 


Assuming that the flow field incident to the orifice consists of a 
steady-state part and an oscillating part, then we can write 

- f * + e,t‘) 

f) . ^ e, f) 


Substituting Eqn. (A-4) into Eqns. (A-1), (A-2), and (A-3) and 
subtracting out the steady-state terms yields 
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To be precise_ii* and v* should be replaced with the expressions 
y* + and V* + v*^ but this refinement will be ignored. Equations 

(7) , (8) and (9) of Section 2 follow immediately upon substituting 
into Eqns. (A-5), (A-6) and (A-7), the non-dimensional variables 
r, t, u, p, p from Eqn. (6) and the characteristic values of P*, 

P*» V*, (D*+L*) and described in Section 2.1. 

Solution to Eqn. (116) 

Introducing the transformations F(t;8) = /2" G(t;6), 6 = /2g 
into Eqn. (95) results in a simplified equation 

5 4t- + <3^ = Cos(t') (A- 8) 


Since we are interested in the behavior of Eqn. (A-8) in the neigh- 
borhood of t = ir/2 (where it is singular - see discussion following 
Eqn. (97), we introduce the transformations 

= QC'Cj ^) , r^Vz-t (A-9) 


into Eqn. (A-8) resulting in 

-S ^ + 6* = sin (t) 

dV 


(A-IO) 
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We are interested in the behavior of Eqn. (A-IO) near t = 0. An- 
ticipating a boundary- layer type behavior near x = 0, the solution 
to Eqn. (A-10) is divided into two regions, one called the "outer" 
region where 

CA-ll, 


and the other the **inner** region where 


I dt 



(A-12) 


An approximate solution to Eqn. (A-11) for 6=0 and t small is 

G»(Tj^=o) ^ Vsin('t') ^ 


In the terminology of singular perturbation theory, Eqn. (A-13) is 
called the outer solution. Note that dG/dx is infinite (singular) 
at X = 0 (t = ir/2) . To remove the singularity, the inner variable 
5 and inner function f(5;<S) are introduced 




where a and b are constants that are determined below. 


The idea here is to transform both the dependent and independent 
variables of Eqn. (A-10) in such a way that both terms on the LHS are 
of equal importance. The transformation is valid only within a small 
region near t = 0 (or t = 7r/2). The boundary condition governing 
the transformed inner solution are determined by proper matching for 
large values of 5 to the outer solution given by Eqn. (A-13). To 
be precise, proper matching occurs when the inner behavior of the 
outer solutions matches asymptotically with the outer behavior of 
the inner solution. ^To demonstrate this, Eqn. (A-13) is rewritten 
with T replaced by 6"^^, 




gV. ^ (A-IS) 


38 



II 


The transformation defined by Eqn. (A-14) must asymptotically (for 
C large) match Eqn. (A-15) . Substituting Eqn. (A-14) into Eqn. 
(A-10) and expanding for small S, t fixed, yields 




(A-16) 


A meaningful solution to Eqn. (A-16) follows upon setting 

I + b -o =2b = a 

yielding 

a = Va , b = '/a 

Thus 

Q (t ^ -f ( ^ j i 

Substituting Eqn. (A-18) into Eqn. (A-16) yields 

+ 4-r/3(i4’).o(r/o 

To solve Eqn. (A-20), expand f as follows 

where the function fgC?) and f^CC) must have the following asymptotic 
behavior to properly match the inner and outer solutions. 

(A- 22) 

4->oo 

and 

•P.fi') = - -jk (*-21) 


(A-17) 


(A-18) 


(A-19) 


(A-20) 
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Substituting Eqn. (A-21) into Eqn. (A-20) and collecting the 
coefficients of 6® and 6' yields 
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Orifice Outlet Velocity 



FIGURE 1. INGARD AND ISING EXPERIMENTAL ARRANGEMENT 
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FIGURE 3. THURSTON'S ET. AL. MEASUREMENTS OF ORIFICE IMPEDANCE 
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FIGURE 4A. COORDINATE SYSTEM USED IN ORIFICE MODEL 
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FIGURE 5. COMPARISON BETWEEN PREDICTED 
AND MEASURED ORIFICE IMPEDANCE 


46 








50 




